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ABSTRACT 



The broad 30 yum feature in carbon stars is commonly attributed to MgS dust particles. However, reproducing the 30 yum feature with 
homogeneous MgS grains would require much more sulfur relative to the solar abundance. Direct gas-phase condensation of MgS 
occurs at a low efficiency. Precipitation of MgS on SiC precursor grains provides a more efficient formation mechanism, such that the 
assumption of homogeneous MgS grains may not be correct. Using a Monte Carlo-based radiative transfer code, we aim to model 
the 30 yum feature of the extreme carbon star LL Peg with MgS dust particles. We find that for LL Peg this modeling is insensitive to 
the unknown MgS optical properties at A < 10 /jm. When MgS is allowed to be in thermal contact with amorphous carbon and SiC, 
the amount of MgS required to reproduce the strength of 30 yum feature agrees with the solar abundance of sulfur, thereby resolving 
the reported MgS mass problem. We conclude that MgS is a valid candidate to be the carrier of the 30 ^m feature when it is part of a 
composite grain population that has optical properties representative of an ensemble of particle shapes. 
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1. Introduction 

The broad 30 yum feature in the thermal c ontinuum emission 
of car bon stars was identified in 1985 by [Cjoebel & Moseleyl 



1985 ) as due to magnesium sulfide (MgS). Begemann et al. 



19941) presented optical data for MgS and compared them to 



the thermal continuum emission of CW Leo to confirm MgS 
as the likely carrie r of the 30 um feat ure. The Infrared Space 
Observatory (ISO; iKes sler et al. 1996h ob served a large sam- 
ple of carbon stars (e . g. |, Yamamii ra et al.l 1998; Jiang et al. 1 999; 
ISzczerba etaDll999t iHrivnak et".! 120001: IVolk et al.ll2002l) dis- 



playing a diversity in strengt h, shape and width of the feature. 
In an analysis of this data set. iHonv et al.l (l2002h concluded that 
the shape of the 30 pm feature c an be b est reproduced when 
the grains are not perfect spheres. IZhang et al. (2009) modeled 
the 30 pm. feature in the proto-planetary nebula HD 56126 using 
pure MgS grains. Assuming the grains to be irradiated by unat- 
tenuated stellar light, their analysis needed the optical properties 
of MgS at wavelengths A < \Q pm. Unfortunately, such mea- 
surements are lacking, as yet. Assuming relatively high absorp- 
tion efficiencies in this regime, these authors required an amount 
of MgS up to ten times the amount of available atomic sulfur 
to explain the strength of the 30 yum feature. Owing to this mass 
problem, Zhang and collaborators argued against MgS as the car- 
rier of the feature. 

In this letter, we report on the observational evidence for 
composite grains in the outflow of the high mass-loss rate 
asymptotic giant branch (AGB) star LL Peg, also known as 
AFGL 3068. We show that MgS is a viable candidate to explain 
the 30 pm feature independent of the absorbing efficiencies of 



these grains at A < 10 pm. Moreover, we show that if one as- 
sumes thermal contact between the dust speci es in the outflow , 
the mass problem in HD 56126 reported by iZhang et al.l (l2009h 
does not occur in LL Peg. Finally, we discuss the formation of 
composite grains in AGB outflows. 



2. Data 

The spectral energy distribution (SED) of LL Peg was con- 
structed by combining spectra taken with the Short Wavelength 
Spectrometer (SWS) and Long Wavelength Spectrometer (LWS) 
on board ISO. The SWS data were retrieved from the Sloan et al] 
(2003) database; the LWS data from the ISO data archive. The 
LWS data were rescaled to the calibrated SWS data. We cor- 
rected for interstellar re ddening following the extinction law of 
IChiar & Tieleiisl (l2006h . with an extinction correction factor in 
the K-band of Ak = 0.01 mag dPrimmel et al...2003.) . 



3. Modeling the thermal energy distribution 

We combined two numerical codes to repro duce the ISO ob - 
servations of LL Peg. The first is MCMax dMin et al.l l2009l) . 
a dust continu um radiative transfer cod e based on a Monte 
Carlo method (iBiorkman & Woodll2001h . which computes the 
dust temperature self-consistently from the thermal energy bal- 
ance equation. The second is GASTRoNOoM (Decin et al. 20Q^ 
I20T0I) . which calculates the momentum transfer from dust to 
gas in the outflow. The results of these two modeling tools 
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ISO - SWS and LWS 

— no thermal contact [CDE a-C] 

— no thermal contact [DHS a-C] 
thermal contact [DHS a-C] 




A (/urn) 



Fig. 1. SED of LL Peg. The SWS and LWS data are shown in full 
black. The dashed red curve shows a model assuming thermal 
contact between all dust species, whereas the full blue line shows 
a model assuming no thermal contact. The dust composition of 
both models consists of a-C grains in DHS shapes and SiC, MgS 
and Fe grains in CDE shapes. The full green model includes a-C 
grains in CDE shapes and assumes no thermal contact. 
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Fig. 2. Dust temperature profile of the envelope of LL Peg for 
individual species and composite grains: Fe in red, a-C in blue, 
MgS in dashed-black, SiC in dashed-green and composite grains 
in magenta. Note that the red, blue and magenta profiles essen- 
tially coincide throughout the whole envelope. 

were iterated to achieve a sel f-consistent solution for the ob- 
served gas expansion velocity dLombaert et al.ll2"012h . Recently, 
iDe Beck et all ilQldi found Mg = 3.1 x 10"^ Mq yr ' for LL 
Peg, but indicated a possible mass-loss variability. We assumed 
Mg = 3 X 10"^ Mq yr"', which leads to a dust expansion ve- 
locity of 16 kms"'. The stellar effective temperature of LL Peg 
was taken to be Ti, - 2400 K. The adopted temperature does not 
impact the results because the CSE of LL Peg is optically thick 
up to /i ~ 50 jum. We used a luminosity of - 1.1 x 10"* L© 
at a distan ce of dj, - 1300 pc, following the period-luminosity 
relation of lGroenewegen & Whiteloc3 (Il996h . The stellar radius 
is thus = 608 Rq. 

Table [1] lists the dust species used in this study and relevant 
dust properties. The bulk of the dust composition is provided by 
amorphous carbon (a-C). We included a small fraction of metal- 
lic iron (Fe), which is a common dust c omponent in the out- 
flow of carbon stars dLattimer et alJ[T978h . We used silicon car- 
bide (SiC) to model the 11 jum feature (see ISpeck et al.ll2009l 
and references therein for a more detailed discussion on SiC 
as a dust component in carbon-rich envelopes), and MgS to 
model the 30 fim feature (see Sect. |4|i. The optical properties 
of these dust species were derived from reflection measurements 
in the laboratory. References for these measurements are listed 
in Table [H We used three models to rep resent grain shapes 
(iBohren & Huffman! [T98l iMin et al.l 12001 : a continuous dis- 
tribution of ellipsoids (CDE), a distribution of hollow spheres 
with filling factor 0.8 (DHS), and spherical particles (MIE). We 
adopted a single grain size of 0.01 /im for all grain shape models. 

In constructing the optical properties of composite grains 
we summed up the extinction contributions of the separate dust 



Table 1. Chemical formula, specific density ps, condensation 
temperature T^ond, and condensation radius 7?cond of the dust 
species. The final column li sts the referenc e to the optical con - 



st ants: l.lJageret all (Il998l): 2. iH enning & St ognienkol (Il996l) : 



3. lPitman et al.l(l2008h : 4. lBegemann et al.i(il994l) 



Dust species 


Chem. 


Ps 




^cond 


Ref. 




form. 


(g cm" 3) 


(K) 






Amorphous carbon 


a-C 


1.8 


1650 


2.5 


1 


Metallic iron 


Fe 


7.9 


1100 


7.0 


2 


Silicon carbide 


SiC 


3.2 


1400 


4.0 


3 


Magnesium sulfide 


MgS 


3.0 


750 


11.0 


4 



species. This is a valid assumption if absorption and emiss ion oc- 
cur in the Rayleigh limit, i.e. at » 0.01 iJ.m (iMin e t al 200|). 
Moreover, the refractive indices of MgS are significantly higher 
in the 30 yum feature than those of the other dust species included 
in the model and, for the relative abundances found in this study, 
dominate the spectral behavior of the composite grains in this 
wavelength region. This supports our approach. Alternatively, 
one may compute the optical properties of composite grains. 
This would require assuming a composite structure, which may 
add a larger uncertainty to the spectral shape than the assumption 
of summing the individual extinction contributions. 

The 1 1 yum feature is best fitted with SiC dust in a CDE en- 
semble. Using the CDE ensemble for a-C grains, however, sig- 
nificantly increases the emission at /I > 50 fim, as shown by the 
full green curve in Fig.[T] The full blue model, on the other hand, 
adopts the DHS particle shapes for a-C, yielding better results. 
We find a dust-mass-loss rate of Md = (1.7+0.1)x 10"^ Mq yr"', 
comprised of 70% a-C, 5% Fe, 10% SiC and 15% MgS. 



4. The 30 yum feature: resolving the mass problem 

In principle, the result of a model including MgS may critically 
depend on the optical constants assumed for MgS at A < 10 
fim, since the shape of the 30 jum feature is very sensitive to the 
temperature distribution of the MgS grains (jHonv et al. 200^. 
For stars that have a low mass-loss rate this introduces a large 
uncertainty in the spectral behavior of MgS because the heating 
of this species is caused by direct stellar light at precisely these 
wavelengths. However, for a high mass-loss-rate star such as LL 
Peg, the model dependence on the unknown short wavelength 
MgS optical constants disappears. The wind of LL Peg is opti- 
cally thick up to /I ~ 50 yum with the transition from an optically 
thick to an optically thin envelope in the infrared occurring be- 
tween 15 and 35 R*. For individual species the dust temperature 
profiles coincide in the optically thick region of the outflow, see 
Fig.|2l In the optically thin, outer wind the MgS and SiC profiles 
start to deviate from those of a-C and Fe. The 30 yum feature is 
produced by emission from dust particles in the outer envelope, 
which are heated by the infrared radiation emitted from an op- 
tically thick surface at ~ 25 R*. The dust temperature at this 
radial distance is ~ 400 K, such that most of the heating occurs 

at /I > 7 yL(m. 

As indicated by IZhang et all (l2009h . the strength of the 30 
yL/m feature in HD 56126 cannot be reproduced using pure MgS 
grains, unless the MgS abundance far exceeds the amount of 
atomic sulfur that is available from atmospheric abundance es- 
timates. This mass problem is shown in the full blue model in 
Fig.[T] which assumes an amount of Mg S dust in agre ement with 
the solar abundance of sulfur However. IZhang et al.l used chem- 
ically homogeneous grains, which heat and cool independently 
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of grains of a different chemical composition. If all dust species 
are in thermal contact, the temperature distribution (shown by 
the magenta curve in Fig.|2]i is predominantly set by the chem- 
ical component that is both reasonably abundant and is an effi- 
cient absorber For a carbon-rich environment, a-C grains heat 
and cool more efficiently over a broader wavelength range than 
MgS or SiC grains. As shown by the dashed red curve in Fig. [1] 
assuming thermal contact between the dust species, but other- 
wise identical parameters as in the blue model, indeed results in 
a 30 fim feature that is produced well in both strength and shape. 

We quantified our results by calculating the atomic num- 
ber abundance with respect to H2 of sulfur needed for the MgS 
dust mass in our models. We only considered sulfur, since it is 
the least abundant component with a so lar atomic abundance of 
(2.6 + 0.2) X 10"^ dAsplund et al.ll2009h . Two models were cal- 
culated, one assuming thermal contact between dust species (see 
the dashed red curve in Fig. [TJ, and one assuming no thermal 
contact (see the dashed blue curve in Fig.[3]for the continuum- 
divided 30 //m feature). In both models, the dust composition 
and dust-mass-loss rate were adapted such that the 30 /im fea- 
ture is reproduced with a similar equivalent width. The model 
assuming thermal contact yields an atomic number abundance 
of 2.5 X 10"^ for sulfur The model without thermal contact re- 
quires an abundance of 1.2 x 1 Given the potentially vari- 
able mass loss of LL Peg (e.g. iMauron & Huggirisll2006l) . we 
estimate the uncertainty on these values to be a factor of two. 
In the case of thermal contact, the required amount of MgS to 
reproduce the 30 /im feature is within the limits imposed by a 
solar sulfur abundance, independent of the particle shape model. 
Without assuming thermal contact, the required amount of MgS 
would significantly exceed a solar atomic sulfur abundance. 

5. Discussion 

5. 1 . Homogeneous versus composite grains 

The above motivated need for thermal contact between the 
dust species implies that they must be included in some 
kind of heterogeneous c omposite grain structure. Recently, 
IZhukovska & Gaiil (l2008h discussed the possibility of forming 
MgS dust in carbonaceous environments through precipitation 
on SiC precursor grains, which would be a plausible way to 
achieve composite grains. These authors showed that the for- 
mation of MgS is strongly coupled to that of SiC, a process 
in which sulfur is being freed by breaking up SiS molecules. 
This sulfur is stored in H2S, which can react with freely avail- 
able magnesium to condense into MgS. A connection between 
molecular SiS and th e formation of MgS is also supported by 
ISmolders et al.l (l2012h . who found a strong correlation between 
the presence of SiS molecular bands and a 30 /im feature in their 
sample of S-stars. SiC is not formed in these stars, but SiS can 
react directly with Mg to form MgS, albeit less efficiently. This 
scenario agrees with the significantly lower MgS abundance that 
is needed to explain the typically weak 30 fim feature in S-stars. 

Additional s u pport for the above hypothesis was given by 
iLeisenring et al.l (l2008h . who compared galactic carbon stars 
with a sample in the LMC, in which all sources show an SiC fea- 
ture at 1 1 yum, while only half of the sources show emission at 
30 fim. Stronger 30 fim features are found together with weaker 
SiC features, which led Leisenring and collaborators to suggest 
that MgS forms as a coating on top of SiC grains. In galactic 
AGB stars, a-C is expected to form first, followed by SiC. If it is 
energetically beneficial for SiC to form on top of a-C grains in a 
core-mantle structure, or for homogenous a-C and SiC grains to 




A (/y,m) 

Fig. 3. Continuum-divided 30 ftm feature with the S WS data of 
LL Peg shown in full black. The CDE model in full red assumes 
thermal contact and requires a sulfur abundance of about the so- 
lar value. In the model represented by the dashed blue line the 
CDE particles are not in thermal contact. To fit the strength of the 
feature sulfur needs to be ~ 5 times the solar value. A DHS (full 
yellow line) and MIE (dashed green) model represent alternative 
shape distributions of particles in which the dust components are 
in thermal contact. The sulfur abundance in the MIE model has 
been adjusted to match the equivalent width of the observed fea- 
ture. 

stick together in an aggregate structure, it is likely that compos- 
ite gra ins are formed. For core-mantle grains. IZhukovska & Gaill 
(l2008h showed that a resonance effect caused by such a structure 
induces a second peak at ~ 35 /im in the absorption efficiency 
profile of spherical MgS grains. However, an ensemble of non- 
spherical grains will likely not produce such a resonance effect, 
see Sect. 15.21 

MgS formation directly from the gas-phase - without any 
kind of precursor grain of different composition - would pro- 
vide a mechanism to produce chemically homogeneous grains. 
However, several studies have indicated that this process is 
not sufficiently efficient to explain the lar ge amounts of MgS 
needed to produce a stro ng 30 //m feature (iKimura et al.ll2005l: 
IZhukovska & Gaill 120081: ICherchneffI [20111). In their study of 
carbon-rich AGB stars and planetary nebulae. pHonv et aT1(l2002h 
adopted chemically homogeneous grains and derived a temper- 
ature for the MgS particles by fitting the 30 /im feature after 
subtraction of a smooth continuum. For LL Peg, they found 
a continuum (i.e. essentially a-C) temperature of 340 K and a 
MgS temperature of 120 K. They ascribed the difference to the 
chemical homogeneity of the grains. However, the authors as- 
sumed that MgS is formed in an optically thin medium and can 
be characterized by a unique temperature. Our radiative trans- 
fer calculations show that both these assumptions are not valid 
for LL Peg . Moreover, the low temperature of MgS found by 
iHonv et all leads to the MgS mass problem. We conclude that 
a population of chemically homogeneous grains, and hence no 
thermal contact between dust species, cannot be reconciled with 
the spectrum of LL Peg. 

5.2. Particle sliape and size 

The spectral shape of a dust emission feature may strongly 
depend on t he model that is used to describe th e par- 
ticle shapes (iHonv et al.l 120021: iMutschke et"an l2009l) . The 
continuum-divided MgS features for several particle models 
are shown in Fig. [3] Mie particles clearly cannot reproduce 
the spectral profile around 30 urn in LL Peg. Several studies 
(e.g. IHonv etaljl200l iMin et al.lllool have indicated that Mie 
particles in any grain-size distribution cause narrow features and 
certain resonances, which have not been observed in the thermal 
continuum emission of AGB stars. It is generally better to use an 
ensemble of particle shapes, such as given by a CDE or DHS 
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Fig. 4. Continuum-divided ISO spectra of multiple carbon stars 
in the region of the 30 fim feature, with LL Peg as a reference 
in black. The top panel shows three typical AGB stars with rela- 
tively low-density envelopes; the bottom panel shows three typi- 
cal AGB stars with rela tively high density-envelopes (data taken 
from lHonv et"al]|2002l) . 

(see iMutschke et al.l |2009| for an overview of these and more 
advanced particle shape models), because they more accurately 
represent features caused by a collection of irregularly shaped 
dust particles. Note that for both the CDE and the DHS model 
the extinction properties are independent of our single grain-size 
assumption as long as the Rayleigh limit holds. The assumption 
of homogeneous or composite grains also impacts the shape of 
the MgS feature, because the dust temperature profile is affected. 
If one assumes homogeneous grains (i.e. no thermal contact) and 
increases the MgS-dust mass such that the strength of the feature 
is approximated, one finds that the 30 jim feature shape is rather 
poorly reproduced by a CDE model (see the dashed blue curve 
in Fig. [3]). However, this poor agreement cannot be used as an 
argument for the need of thermal contact because CDE is quite 
a simple grain shape model. A comparative study between par- 
ticle shape models bv lMutschke et al.l (l2009l) has shown that the 
CDE model often results in emission bands too much enhanced 
at red wavelengths when compared to dust optical properties 
measured in the laboratory. The good fit of our MgS CDE model 
given by the red curve and the observed 30 yum feature in LL Peg 
consequently may indicate some additional effect that enhances 
the red wing of the observed band profile. Condensation experi- 
ments bv lKimura et"aD(l2005h have shown that MgS dust formed 
as network-like structures can show an enhanced red part of the 
MgS band profile. Spectra of dust grains embedded in a matrix 
also show such an enhancemen t compared to those of free par- 
ticles (e.g. iTamanai et al.ll20()6l) . which may point to mixing or 
embedding of the MgS dust in inhomogeneous grains. 

5.3. Diversity of the 30 /im feature shape in AGB outflows 

LL Peg is a highly evolved AGB star and exhibits one of the 
reddest thermal emission continua of all carbon stars. The out- 
flow of this source is optically thick even in the 30 /im fea- 
ture. As indicated in the bottom panel of Fig. H] AGB stars 
with a high-density envelope show an almost uniform shape 
of the 30 fim feature. Sources with a lower density envelope, 
such as those shown in the top panel of Fig. |4] display a wide 
diversity i n shape. As a first p ossible explanation of this phe- 
nomenon, iKimura et al.l (l2005h suggested that the shape of the 
30 fim feature may depend on the formation history of MgS. 
If the grains have formed through gas-phase condensation, the 
resulting emission feature is expected to be more pronounced 



at /l ~ 32 yum than at longer wavelengths. If MgS grains are 
produced via a gas-solid reaction, the red wing of the feature 
at /I > 35 i-im becomes more intense relative to the strength 
at /i ~ 32 fim. As shown in Fig. ID the 30 /im feature in the 
high-density envelopes has a more pronounced red wing than 
those in the lower density envelopes. One may speculate that 
in low-density sources, where the 30 /im feature is optically 
thin throughout the entire envelope, the feature may be com- 
posed of gas-phase condensates in the inner part of the wind 
and grains formed by gas-solid reactions throughout or in the 
outer part of the wind. In high-density sources, one only ob- 
serves the contributions from the latter formation channel. This 
scenario c ould be consistent wit h the MgS-SiC correlation in S- 
type stars (ISmolders et al.ll20T2l) . where the lack of SiC particles 
indicates a gas-phase condensation channel, albeit not a very ef- 
fective one. A second explanation might be that the low opac- 
ities in low-density envelopes favor the heating of another dust 
species with an absorption coefficient profile peaking at about 
30 fim. Alternative carriers sugge sted in the literatur e include 
hydrogenated amorphous carbons dGrishko et al.ll200"TI) . 

6. Conclusions 

We have modeled the ISO spectrum of the high-density carbon 
star LL Peg with a dust composition consisting of a-C, Fe, SiC, 
and MgS. The (high) density and temperature structure in the 
envelope of this source allow one to model the 30 fim feature 
with MgS dust particles, independent of the unknown MgS opti- 
cal properties at A < 10 //m. We showed that MgS is a viable 
candidate to be the carrier of the 30 fim feature. An ensem- 
ble of particle shapes works significantly better than spherical 
grains to explain the shape of the feature. Thermal contact be- 
tween the dust species is required to ensure that the amount of 
MgS dust in the envelope of LL Peg does not exceed the solar 
abundanc e of sulfur, thereby avoiding the mass problem as re- 
ported bv lZhanget aH (l2009h for the post-AGB star HD 56126. 
Achieving thermal contact between all dust species is possible 
if these species form in some kind of heterogeneous composite 
grain structure. 
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